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Abstract 

We analyze the photon statistics of a weakly driven cavity qnantnm elec¬ 
trodynamics system and discnss the effects of photon blockade and photon- 
induced tunneling by effectively utilizing instead of avoiding the center-of- 
mass motion of a two-level atom trapped in the cavity. With the resonant 
interaction between atom, photon and phonon, it is shown that the bunching 
and anti-bunching of photons can occur with properly driving frequency. Our 
study shows the influence of the imperfect cooling of atom on the blockade 
and provides an attempt to take advantage of the center-of-mass motion. 

Keywords: Cavity quantum electrodynamics, photon statistics, 
center-of-mass motion 


1. Introduction 


The nonlinear interaction represents a hot research topic in recent years. 
It can lead to a variety of intriguing nonlinear optical phenomena which are 
ubiquitous in quantum optics, such as the vortex formation, self-focusing, 
soliton propagation, optical multi-stability and so on 0,3- The nonlinear 
interaction between photons, light or under the influence of Kerr medium 
3 3 has also been widely used to generate nonclassical field states 0 3 
and in the observation of strict quantum effects. Cavity quantum electro¬ 
dynamics (CQED), as the light and matter meeting interface, is the most 
straightforward way to produce the quantum nonlinear dynamics BB- A 
lot of strict quantum effects such as optical solitons 0, quantum phase tran¬ 


sitions [10, [11|, quantum squeezing [12[ and optical switching with single 
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photon [13[ have been demonstrated in CQED systems based on the strong 
optical nonlinearity. 

As a typical nonlinear quantum optical effect, photon blockade shows that 
the system ‘blocks’ the absorption of a second photon with the same energy. 
The typical feature of this effect is the photon anti-bunching which is signaled 
by a rise of with r increasing from 0 to larger values while < 

9'">(r) as discussed in detail in Ref [1^ . Recently, quantum blockade schemes 
have been reported in various physical setups such as sup erconducting circuit 
, cavity arrays [ll, l3^19|, quantum dots Mm, atomic systems [3, 


conhned cavity polaritons 
What h 


IS more 


I 22 L |23| . quantum optomechanical setup s |24l -l27 
28[, systems with ultrastrong coupling |29| and so on 
, the strong bunching effect, and its special signature called photon-induced 
tunnelling, which shows the absorption of the hrst photon can enhance the 

3- As to the 


absorption of the subsequent photons are also been studied 
photon blockade in a cavity, the fundamental mechanism is the anharmonic 
energy-level structure of the light hel^when the photon-photon interaction 

or induced by the 


is induced by the nonlinear medium 
interaction between the trapped multi-level atom and the multi-mode cavity 
^1, or is observed in the detuned Jaynes-Cummings (J-C) model [^, 42 . 

In this paper, we study the photon statistics in the detuned J-C model 
and the effect of the atomic center-of-mass (COM) motion on the photon 
statistics. Instead of the multi-level atom and the multi-mode cavity, we re¬ 
quire that a single two-level atom trapped in a single mode cavity driven by 
a weak laser and the atom oscillates at its origin with trap frequency. Thus, 
the photon-photon interaction is induced by the participation of phonons. It 
is found that the photon blockade is generated under the resonant interaction 
among atom, photon and phonon, and both the bunching and anti-bunching 
of photons strongly depend on the atomic motion. In particular, by com¬ 
paring the cases with and without atomic center-of-mass motion, we find 
that the atomic center-of-mass motion can enrich the phenomenon of photon 
statistics. The advantages: on one hand, are to reveal the influence of the 
atomic center-of-mass motion on the blockade due to the imperfect cooling 
and on the other hand, are to provide an attempt to taking advantage of the 
center-of-mass motion of the trapped atom. This paper is organized as fol¬ 
lows. In Sec 2, we analyze the equal-time correlation function of cavity mode 
and demonstrate the photon statistics of J-C model; In Sec 3 we introduce 
the atomic center-of-mass motion, and investigate the equal-time correlation 
two-time correlation function of photons; Finally, the conclusion is given in 
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Sec 4. 


2. Equal-time correlation in the detuned J-C model 


A two-level system (ion or an atom) couples to the cavity with frequency 
Ua which is driven by an external optical held. The frequency of atomic tran¬ 
sition from ground state \g) to excited state |e) with linewidth 7 is denoted 
by Ue- The related Hamiltonian under the rotating wave approximation can 
be written as 

H = Aa^a -|- 6a^a~ + g{a^a~ + aa~^) -|- (a^ -|- a) , (1) 


where A = Ua — ojl is the laser detuning from the cavity mode and 6 = 
ul is the laser detuning from the atom, in the n-excitation subspace, 


OJ. 


the eigenstates and the eigenvalues can be easily get l44 


n+) 



( 2 ) 


n—) 



( 3 ) 


with 6 = 6 — A the detuning in frequency between the two-level transition 

system and cavity mode and A'=\JAng'^ + 5 the energy difference between 

manifold levels. The eigenvalues are = {n — |)A + A ^{n ^ 1). 
We can easily hnd that the energy splitting has a nonlinear dependence on 
n (photon number) and gives rise to the nonlinear optics which contains the 
photon blockade effect and tunneling phenomenon. We will show this in the 
following part by a strict treatment. 

Fig. 1(a) shows the energy diagram of the system. The energy difference 
value A' which is a not a constant, this anharmonic spacing level can affect 
the photon statistical distribution, such as the photon blockade and photon- 
induced tunnelling. To understand the photon statistical situation, we will 
focus on the photon correlation that characterizes the nonclassical photon 
statistics in the system. Here, we study the equal-time (namely zero-time- 
delay) second-order photon-photon correlation function j44l |: 
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Figure 1: (color online).(a) Diagram for the eigensystem of the Hamiltonian of the coupled 
cavity-atom system. The anharmonic energy difference between consecutive manifolds is 

A' = \JAng"^ + S which is not constant, (b) The energy diagram with the anharmonic 
spacing energy under the strong coupling conditions (parameters are chosen as follows). 
The level shift is caused by the strong coupling. The blue lines indicate the frequency 
transitions of single photon process for the |0) —^ |1±) and red line indicate the two- 
photon process |0) —^ |2±) . (c) ,(d) Mean photon numbers of the cavity and the equal-time 
second-order function vs. the cavity-laser detuning A for the J-C model and atomic center- 
of-mass motion model. The red and purple curves are approximate analytical solutions of 
the mean photon number of the J-C model and COM model respectively, while the blue 
and green dashed curves stand for numerical results based on the solution of the quantum 
master equation Eq. (13) and Eq. (32). Further explanation of the difference about the 
two model can be found in the next section. Here all parameters are dimensionless. The 
parameters are 7 /k = 1, F/k = 0.1, gj k = 50, J = A -|- J, H/k = 0.1, J/k = 50. 
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where n = {a)a) is the intra-cavity photon number of the cavity mode, Pn 
represents the probability with n photons. In Eq. (jl]) the operator is evalu¬ 
ated at the same time. When the photon anti-bunching occurs, the second- 
order correlation function should fulhlls the inequality < 1 and the 

limit —)■ 0 corresponds to the perfect photon blockade in which two 

photons never occupy the cavity at the same time. On the contrary, when 
> 1, it means photons inside the cavity enhance the resonantly enter¬ 


ing probability of subsequent photons j4a, 146 


In order to give an intuitive picture, we take an analytical (but approxi¬ 
mate) method to calculate the second-order correlation function by employ¬ 
ing the wave function amplitude approach. Considering the effects of the 
leakage of the cavity k, the spontaneous emission 7 of the atom, we phe¬ 
nomenologically add the relevant damping contributions to Eq. ([1]). Thus 
the Hamiltonian can be rewritten as H — i{Ka^a + a~). Since we consider 
the weak driving limit, only few photons can be excited in the cavity. So one 
can assume that the state of the composite system can be given by 


48 


I'J') — |0, S') + Hig |1, gf)-f Hoe |0, e) 

-\-Aie |1, e) -1- A 2 g |2, g). 


(5) 


Since the dynamics of the system is governed by Schrodinger equation, using 
the Hamiltonian and the state |4/) , we can arrive at the equations about the 


amplitudes in Eq. (jS]) as follows. 

Aig = —(k -|- i/A)Aig — iVtAog — igAoe — \/ 2 iQA 2 g, (6a) 

Aoi =-{'f+ i6)Aoi-igAio-inAii, (6b) 

A2g = —2{k + i/A)A2g - V 2 igAif,, (6c) 

Ale = —{1^ + 7 + -l- i/S)Aie — V2igA2g — iVLAie. (6d) 


Solving Eq. (6) will reveal all the physics. Let the initial state of the system 
be lO,^)- Considering the weak limit of the driving field again, we can get 
Aiig —)■ 1, and the Eq. (6) are closed. Thus, Eq. (6) can be easily solved. In 
the following, we will only consider the question in the steady-state case. In 
addition, the steady-state solution of Eq. (6) can be analytically obtained, 
but the concrete form are quite cumbersome, so we further neglect the high- 
order terms in H (weak driving) and obtain 



zr2(7 -f i5) 

^2 + A 


(7) 
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T 9^ 

9= + A’ 

(8) 

n^[g^ - {-f + tdf - A] 

V2{g^ + A){g^ + {k + tAy + AY 

(9) 

ig‘^Vt^{'y + ti + id + iA) 

{g^ + A){g^ + {k + lAf + AY 

(10) 


with A = {'y+i6){K+iA). Thus Eq. (jl]) can be rewritten as ’ 

with Pi = \Aig\ ,p 2 = \A 2 g\ ■ In the weak-driving case, we can easily get 
Pi 3> P 2 , then the second-order correlation function 




+ y + - {'y + i5Y - A 


I 7 -f + {k + iAy + A 
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( 11 ) 


I 7 i6\^ 

g^ + A^ 


( 12 ) 


In order to show the validity of the above results, we use the Markovian 
master equation for the model, that is. 


p=-i[H,p\ +KL[a]p + 'yL[a ]p (13) 

where H is the Hamiltonian given by ([T]), p is the density operator of the 
whole composite system, and L[d]p = 2dpd^ — (Adp — pd^d, {d = d,d~) is the 
dissipator. In addition, we don’t consider the thermal photons for simplicity 
( 4 ^ . Since the steady-state solution is needed for our purpose, we will directly 
employ a numerical way to solving Eq. fITT]) for the steady state ps ji^. So 
the second-order correlation function can be directly obtained by g^^^(O) = 

jTr(pslfa)p ’ mean photon number is given by n = Tr{psa^a). 

In what follows, we will employ both the analytical method given by Eq. 
(ITT]) and the numerical way as Eq. f[T5]) to study the properties of the equal¬ 
time correlation function As is shown in Fig. 1(d), we plot p*'^^(0) as 

a function of the cavity-laser detuning A, a different story takes place with 
resonant condition when the model under the strong coupling region; the 
detuning A/k = —50(25) corresponds to the red lines in Fig. 1(b), which 
means the transition |0) —)■ |2±) and p^^^(O) S> 1, indicating that the photon 
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bunching for the cavity. In addition, when A/k = —50,it also indicates a 
quasi-Dark state ( d) (x g |0, 10, e)), which exhibit the strong bunching 


behavior. At the point of A/k = —25(1 ± \/5), 1, indicating the 

complete photon blockade due to the suppressed two-photon process, corre¬ 
sponds to the blue lines in Fig. 1(b). When A = 0, the system does not 
exhibit the strong bunching behavior as the bare cavity j^, that means the 
absorption of the hrst photon can not enhance the absorption of the sub¬ 
sequent photons and the expected photon-induced tunnelling phenomenon 
disappeared owing the cavity-atom off-resonance. 


3. Atomic center-of-mass motion on photon statistics 
3.1. The model and equal-time correlation 


In this section, we consider the effect of atomic center-of-mass motion 
on the photon statistics. As shown in Fig. 2(a), a two-level system which 
is conhned by a harmonic potential with the trap frequency u. In this con- 
hguration, the two-level system and cavity are same as before. We suppose 
that the model is a one-dimensional model and restrict the motion of the 
atom along the x axis. So the Hamiltonian of this system under the dipole 


approximation reads [50 


•^9 9^9 

I _ p mv X 

H = Uaa'a +Uea^a +- -^--— 

2m 2 


+ g{kx){a+a + aV") + H , 


(14) 


where g is the position-dependent coupling coefficient, H and ojl are related 
to the power and frequency of the driving laser, respectively, k denotes the 
wave number of the held. Note that the third and forth terms represent the 
kinetic energy and the harmonic potential respectively, and m represents the 
mass of the atom. The strength of the interaction of the two-level system 
with the single mode of cavity is characterized by the coupling operator which 
is given by |5]| 

g{kx) = gcos{kx cos (j) + if)., (15) 

where 0 is the angle between wave vector and axis of the motion, ip accounts 
for the displacement of the trap center with respect to the origin. For con¬ 
venience we set 0 = 0, p = I , and we also introduce the annihilation and 
creation operators b and of a quantum of vibrational energy, therefore. 
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Figure 2: (color online), (a).The center-of-mass motion of a trapped atom under the as¬ 
sumption that the atom is confined within a high-finesse cavity which is pumped by a 
laser, (b) Diagram of the eigensystem of the Hamiltonian under the resonance condition 
and level diagram is limited in the subspace spanned by the zero-, one-, and two-photon 
states. States are labeled by |na,n&,*) denoting the photon numbers of cavity mode , 
phonon number of mechanical mode and the level of atom. (c),(d): The mean photon 
numbers of the cavity and the equal-time second-order function (0) vs. the trap fre¬ 
quency V respectively. The red curves are approximately analytical solution of Eq. (30) 
and Eq. (31) while the purple curves are numerical results based on the solution of the 
quantum master equation Eq. (32). The parameters are chosen the same as Fig. 1. 



























the position and canonically-conjugated momentum of the atom are given 
by X = + &^), p = —h). In addition, the third and forth 

terms of Hamiltonian flT^ can be rewritten as hv{b^b -l- |). Omitting the 
constants, Hamiltonian flT^ becomes (we set h = 1 hereafter) 

Hq = Uatt^a + uJeO'~^cr~ + vb^b + g{a^a~ + aa~^) (b^ + b^ + 

Q , (16) 


with g = g^y 2 ^- frame rotated at the laser frequency ojl, we obtain 

H' = AaPa + 5a^a~ + vb'^b + g{aP(j~ + aa^) (pp+ b') 

-|-H -|- cPj , (1^^) 

Now, we focus on the three-mode resonant interaction in which the cavity 
and atom exchange a photon by absorbing or emitting a phonon in the me¬ 
chanical mode. We set A = 5 -|- 1 ; and assume 1 In order to study the 

dynamics of the system, we would like to turn to a rotation framework subject 
to the transformation u{t) = exp[—iRt] with R = S{a^a~ — b'^b). With the 
rotating-wave approximation at large detuning 6, the effective Hamiltonian 
can be given by 


hfeff = (5 + v){a^a + b^b) + g{a^ba -|- ab^a^) 

-|- H (a^ -|- a) . (18) 

The second term in the hrst line of Eq. flT8|) describes the effective nonlin¬ 
ear coupling proportional to coupling strength g which describes the coher¬ 
ent photon-phonon exchange between the cavity mode and mechanical mode 
mediated by atomic absorption or emission of a photon. The effective Hamil¬ 
tonian can be diagonalized in the absence of a driving jd^. The eigenstates 
distinguished by different numbers of photons and phonons can be expressed 
as 


|0) = |0,0,^), 

(19) 

l±) = ^(|l,0,^?)±|0,l,e). 

(20) 

2±) = ^(|2,0,^?)±|l,l,e), 

(21) 
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where |na,nf„i) {i = e,g) are the preferential basis with na,nb denoting the 
photon number and phonon number, respectively, and |i) (l^f), |e)) represents 
the ground and the excited states of the atom. The low-energy level diagram 
is sketched in Fig. 2 (b), from which one can see that the states |1, 0, gf) and 
|0, l,e) are superposed to form two eigenstates |1±) splitted by 2g, and the 
states are |2, 0, g) and |1,1, e) are superposed to form another two eigenstates 
|2j-) splitted by 2\/2g. Intuitively, one can hnd, from the level diagram, that 
the resonant absorption of a photon with frequency ujai:g to reach the state 
|1±) ‘block’ absorb the second photon with the same frequency because of 
the detuning from the other levels 0,HS0- 

Using the same approach, we can obtain the following results for the 
atomic center-of-mass motion system. So we can also assume that the state 
of the composite system is 


\^)com = dloo|0,0,r/) + klio(|l,0,r/) + kloi|0,l,e) (22) 

+^11 If, 1, e) + ^20 |2, 0, gf) (23) 

The equations about the amplitudes in Eq. fl2^ as follows. 

Aiq = —[k i/S.)AiQ — iQAoQ — igA^i — \/2inA2o, (24a) 

Aqi = — (7 -|- T -|- zA)74oi — igAiQ — iVtAn^ (24b) 

A 20 = —2(k iA.)A 2 Q — ^/2igAii, (24c) 

All — —(K + yTT-l- 2i/\)Aii — \/2igA2Q — iVtAn^ (24d) 


Solving Eq. (24) will reveal all the physics. Let the initial state of the 
system be |0,0,5'). Considering the weak limit of the driving held again, we 
can get Aqq —)■ 1, and obtain the following coupled set of equations for the 
steady state: 


(7 T zA)yloi 


then we can obtain 


^10 — 


■ ^17^01 ~ \/2iVtA2Q — 0, 

(25a) 

igAio - iVtAii = 0, 

(25b) 

- \/2igAii = 0, 

(25c) 

“ \/2i^yl2o — = 0, 

(25d) 

iQj 

2 -|- 7k ’ 

(26) 
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Figure 3: (color online).We plot the log^Q (0)] as a function of v and 7 for (a) F/k = 0.1 
and for (b) F/k = 1. The other parameters are (5 /k = —50, H/k = 0.1. For (c), we plot 
the logiQ 9 .S a function of v and 6, we take j/k = 1, F/k = 0.1, A = 5 -|- :/,other 

parameters are the same as (a) and (b). 
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Au = 
A 20 = 


g'^ + jK,' 

igkVL^ 


{9^ + 7 /^)( 5 '^ + + 'Jk) ’ 

__ 

'/2{g‘^ + ^k){g^ + + 7 ^) ’ 


where 7 = 7+r+i(5+r^), k = K+i{6+v). Thus withpi = |^iof ,P 2 
Pi ^ P 2 , the second-order correlation function becomes 


(27) 

(28) 
(29) 


A l2 

Ao 


— 1^20 


9f(0) 


2p2 


g^{g^ +^*k*){g^ + ^R) 


Pi {g‘^ + K*'^ +'y*K*){g'^ + k'^ +'yK)\'j\‘^^ 
and the mean photon number is 

{g"^ +'y*K*){g‘^ +'jk)' 

Meanwhile the master equation is 

P = -i[Hes, p] + KL[a]p 'yL[a~]p + TL[b]p, 


(30) 


(31) 


(32) 


where i^eff is the effective Hamiltonian given by flTSll . and L[d]p = 2dpd) — 
(Adp — pd^d, {d = d, b, d”) is the dissipator. 

In what follows, we will employ both the analytical method given by 
Eq. fl30|) and the numerical way as Eq. fl3^ to study the properties of 
the equal-time correlation function ^'•^^^(O). As is shown in Fig. 2 (c), we 
plot 5'*^^^(0) as a function of the trap frequency z/. We hnd that the photon 
statistic properties can be controlled by tuning the trap frequency u. Here, 
the red curves are plotted using the analytical expression Eq. fl30|) while 
the grape curves are based on the numerical solution of Eq. (I32|) . One can 
see that sX v/k = 100(A = 0), (y'*'^^(0) 3> 1 in Fig. 2(d), where the photon 
satisfy the super-Poissonian distribution. In other words, it indicates that at 
this point photons occur the bunching behavior. This can be explained that 
the system presents destructive interference that suppresses the population 
1,0, (yf). At this point, the system is driven into a dark state \dark) cx 


m 


g | 0 , 0 , 5 f) — H |0,1 , e) which is similar to the electromagnetically induced 


transparency In the dark sate, |0, l,e) remains populated, allowing 


transitions to |l,l,e) which is strongly coupled to |2,0,(7). The probability 


12 










of two photons inside the cavity is resonantly enhanced at this point which 
corresponds to the peak in the correlation function of Moreover, 

we see that at vjn = 100, the strong bunching regime is accompanied by a 
suppression of the photon number, the photon-induced tunneling occurs due 
to the probability of single-photon emission decreases, while the probability 
for photon pair generation increases. This behavior can also be occurred 
when vjK = 100 ± 25\/2, because of the transition |0) |2_|_) which is a 

double photon resonance process. At the frequency z//k = 50 and 150, the 
detuning A/k = —50 and 50 correspondingly, the second-order correlation 
function 5 '*'^^( 0 ) 1 (i.e., two dips), which shows the anti-bunching and sub- 

Poissonian distribution. This is consistent with our intuitive analysis given 
above. That is, the frequency of the driving is resonant with the transition 
between the state |0) and |±1), which ‘blocks’ the absorption of the second 
photon due to the detuning. 

So far, we have demonstrated the quantum properties with the analytical 
and numerical methods both in the J-C model and atomic center-of-mass 
motion model. By comparing the two curves in Fig. 1(d) , we hud that 
by introduced phonon degree of freedom, the photon statistics changes as a 
function of A, and photon blockade effect occurs at the different A compared 
with the J-C model which can also be demonstrated by the analytical ex¬ 
pressions of second-order correlation function in Eq. (11) and Eq. (30). We 
can also hnd that the trap frequency v leads to an obvious displacement of 
the maximal peak, whilst it leads to a new bunching peak. 

To hnd the effects of the decay and the detuning on the photon blockade, 
we plot the logarithm of the correlation function 5 '*'^^( 0 ) as a function of 7 
and z/ in Fig. 3. Fig. 3 gives a boundary for different photon distributions 
such as super-Poissonian and sub-Poissonian distributions. The two regions 
are distinguished by y/zt ~ 10 in both Fig. 3(a) and 3(b). These two hgures 
show how the second-order correlation function varies with the decay rates. 
In particular, one can hnd that the maximum of (7*'^^(0) occurs at the top of 
the hgures (y/zt = 15), but the minimum of 5'*'^^(0) (anti-bunching) is just 
below y/zt ~ 10. Between them, there also exists a bunching peak. In Fig. 
3 (c), 5'*'^^(0) between bunching and anti-bunching at the frequency v and 
detuning J is illustrated. It is apparent that the photon blockade happens in 
the region |A/zt;| which can be obtained from Eq. (30) by setting the 
numerator of 5 '*'^^( 0 ) vanishes for small dissipations. 
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Figure 4: (color online).The evolution of finite time delay correlation function 
of the atomic center-of-mass motion cavity mode. The red line indicate the analytical 
expression and the blue dash line denotes the numerical simulation. The parameters are 
chosen as g/n = 20 , Sjn = —70, V,/k = 4, j/k = 4, T/k = 0.1, vjn = 50 . 


3.2. Delayed coincidence correlation 

In addition to the second-order correlation function for equal time dis¬ 
cussed in the previous section, some quantum signatures can also be man¬ 
ifested from the photon intensity correlations with the nonzero delay. In 
this part, we give a discussion of the evolution of the two-time correlation 
function, which is dehned for stationary state by 56 


(2)r ^ (a^(0)at(r)a(r)a(0)) 


(ala)" 


(33) 


Eq. (33) can be rewritten, based on the correlations of classical light intensity 
/, as With the Schwarz inequality for random variables: 

{hhf < (/?> gi) . (34) 


which means that 


{I{T)I{0)f < (/(r)=> (7(0)"). 


(36) 


The inequality is saturated if r = 0. Thus one can further End that [57 1 

<^(2)(o). (36) 


Similar to the classical inequality, > 1 at zero delay, and violation of 

the inequality at hnite delay is an unique signature of quantum system. Next 
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we will calculate the delayed second-order correlation function for the cavity 
mode. The two-time correlation functions can be understood in terms of the 
hve-level model discussed above. For this purpose, we may use the standard 
method to calculate the two-time correlation function with a conditional state 
as the initial condition 58|. The unnormalized state after the annihilation of 


a photon in the cavity is a |T) = AiodO, 0, g) + An |0,1, e) -f \f 2 A 2 Q |1, 0, g). 
By solving the Eq. (ll)(a)-(d) with this state as the initial condition, one 
can obtain the time-delayed correlation function as 



l^io(r)r 


(37) 


The positive derivative of the time-delayed correlation function with respect 
to variable r represents a nonclassical behavior of the system. Additionally, 
it is important for the experimental observation of the anti-bunching, since 
the precise measurement of the equal-time correlation function may be chal¬ 
lenging . Both the analytical and numerical results of the two types cavity 
mode are plotted in Fig. 4. We can observe that according to the system 
nonlinearity the time-dependent second-order correlation function shows the 
quantum signature which violates the inequality in Eq. fl36|l because 
going beyond their initial value at hnite time delay. We note that the be¬ 
havior of the second-order correlation function in this case is similar to the 
single-photon diode in the semiconductor microcavities coupled via x(2) non- 
linearities j^. This can be attributed to the fact the strong nonlinearity of 
our system for both the zero and the hnite time delay correlation functions. 


4. Conclusion 

We have studied the photon blockade effect in CQED system weakly 
driven by a monochromatic laser held, in which an atom is trapped inside. 
The spirit of this scheme is to ehectively take advantage of the atomic center- 
of-mass motion instead of avoiding it. By the numerical and the approximate 
analytical expressions of the one-time and two-time second-order correlation 
function for the cavity photons, we have identihed several diherent processes 
that can lead to photon-induced tunnelling and photon blockade, our study 
provides another way on photon control using a atomic center-of-mass motion 
system. In addition, the trap-frequency dependent blockade ehect directly 
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shows the influence of imperfect cooling, whilst it may induce the other po¬ 
tential applications of atomic center of mass motion in quantum information 
processing. 

Finally, we would like to say that trapped ions or atoms are well suited 
for this purpose as the quantum techno logy for controlling their degrees of 
freedom that is already well developed j^, ^|. The coupling to the cavity 
has been successfully employed for implementation of a photon-mediated 
entanglement distribution j^-^| and the cooling the motion of an atom 
which is trapped by a harmonic trap j^, IgsI-Ig^ . The combination with the 
coupling to the held of high-finesse resonators can bring novel perspective for 
the manipulation of the atomic motion. 
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